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PREPRINT 

DEPENDENCE  OF  MIXED-LAYER  ENTRAINMENT  ON  SHEAR  STRESS  AND  VELOCITY  JUMP 


ABSTRACT 


From  rotating-screen  annulus  experiments  the  entrainment  rate,  wg,  normalized 

by  the  friction  velocity,  u4,  has  been  found  to  be  a  function  of  both  the  overall 

Richardson  number,  RT,  and  the  inverse  Froude  number,  Rv-  The  RT  dependence 

deduced  by  Price  (1979)  and  Thompson  (1979)  satisfactorily  explains  the  present 

-1.4 

data  if  multiplied  by  an  approximate  Ry  ‘  dependence.  The  measurements  indi¬ 
cate  that  R  is  a  variable  that  is  influenced  by  sidewall  friction,  time  after 
v 

onset  of  the  surface  stress,  or  other  factors.  The  greater  w  /u  values  of  the 

© 

Kantha  et  al.  (1977)  type  experiment  over  that  of  the  Kato  6  Phillips  (1969) 
experiment  can  be  explained  by  somewhat  greater  Rv  values  in  the  latter  case. 

A  close  connection  is  now  apparent  between  entrainment  experiments  in  two- 
layer  systems  designed  to  have  only  one  velocity  scale  (the  interfacial  velocity 
jump,  Av) ,  and  the  rotating-screen  annulus  experiments  having  two  velocity  scales 
(uA  and  av) .  The  former  also  have  (at  least)  two  velocity  scales,  the  second  one 
being  associated  with  the  presence  of  turbulence  throughout  one  or  both  of  the 
fluid  layers. 

The  turbulent  layer  is  found  to  be  quite  well  mixed  in  density  only  if 

w  /u^  does  not  exceed  about  0.03,  or  w  /|av|  does  not  exceed  about  0.003.  The 
e  *  e 

present  data  suggest  more  rapid  entrainment  when  temperature  rather  than  salt 
provides  the  density  jump,  as  first  noted  by  Turner  (1968)  in  oscillating  grid 
experiments.  If  this  is  a  Peclet  number  effect,  the  trend  did  not  continue  for 
still  much  greater  Pe  values,  the  data  for  kaolin  (clay)  being  very  compatible 


with  that  for  salt. 
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1.  Introduction 

The  rate  at  which  a  well  mixed  turbulent  layer  entrains  an  adjacent  non-turbulent 
layer  in  geophysically  relevant  situations  is  believed  to  depend  upon  the  surface 
shear  stress, £,  the  depth  of  the  mixed  layer,  h,  and  the  buoyancy  jump  b  across 
the  edge  of  the  layer  (b  =  g  where  g  is  the  gravitational  acceleration,  />0 

a  reference  density  of  one  of  the  layers,  and  J^the  change  in  density  across  the  outer 
edge  of  the  mixed  layer) .  If  a  surface  flux  of  buoyancy  or  heat  is  present,  that 
quantity  may  also  be  important;  however  here  only  the  problem  of  the  neutral  surface 
layer  driven  by  shear  stress  will  be  considered. 

The  first  set  of  experiments  undertaken  to  explore  the  dependence  of  the  growth 
of  a  neutral  mixed  layer  upon  these  factors  was  that  of  Kato  and  Phillips  (1969, 
abbreviated  KP)  using  an  annulus  of  diameter  1.52m  and  gap  Ar  =  0.23m.  The  approx¬ 
imate  relation  they  found  is 

we/u,  =2.5  Rt_1  (1) 

where  w  is  the  entrainment  rate,  or  mixed-layer  growth  rate,  u  is  the  friction 

6  w 

% 

velocity  (|J|/^>o)  ,  and 

Rt  =  bh/u*2  (2) 

is  an  overall  Richardson  number.  Their  actual  data  suggested  a  power  law  some¬ 
what  less  steep  thanR^at  small  RT  values  (RT<50)  and  somewhat  steeper  at  larger 
values  (Rt>100).  In  their  experiments,  using  a  turbulent  layer  of  fresh  water 
(driven  by  a  rotating  screen)  overlying  a  non-turbulent  linearly  stratified  outer  layer 
(abbreviated  SOL)  of  salt  solution,  u*  was  obtained  directly  from  torque  measurements; 
however,  quantitative  velocity  measurements  were  not  made. 

Later,  using  the  same  annulus,  Kantha  et  al.  (1977,  abbreviated  KPA)  explored 
the  two-layer  system  (abbreviated  2LS)  in  which  the  non-turbulent  lower  layer  is 
denser  than  the  turbulent  upper  layer  but  is  not  stratified.  The  entrainment  rate 
was  not  found  to  obey  any  simple  power-law  dependence  upon  Rr;  values  of  w /u4 
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were  about  half  those  found  by  KP  for  the  same  values  of  and  h/ar  (Price,  1979). 

A  significant  reduction  in  w  /u*  with  increased  sidewall  drag  (increased  h/At)  was 
noted . 

Soon  thereafter  further  experiments  with  the  2LS  were  conducted  by  Kantha 
(1978)  using  a  scaled-down  annulus  half  as  large  as  that  of  KP  and  KPA.  His 
w^/u^  versus  Rr  values  closely  resembled  those  of  KPA  but  were  substantially  smaller 
(h/ar  was  generally  larger).  Salinity  profiles  were  also  measured,  and  indicate, 
surprisingly,  that  the  turbulent  layer  often  was  not  very  well  mixed  in  salt  content. 
Although  the  purpose  of  employing  a  smaller  annulus  was  to  determine  if  centrifugal 
effects  were  important,  it  was  not  clear  if  differences  observed  were  associated 
with  such  effects  or  with  generally  larger  sidewall  drag.  Near  the  inner  annulus 
wall,  for  RT>100,  centrifugal  effects  were  apparently  suppressing  the  turbulence, 
judging  from  shadowgraph  observations. 

An  explanation  for  the  discrepancy  between  the  SOL  and  2LS  results  was  provided 
by  Price  (1979)  and  Thompson  (1979).  Price  started  with  the  momentum  budget  for  the 
mixed  layer,  neglecting  centrifugal  effects  and  curvature: 


h(hQ)  =  u*2  -  2cDw'?2h/4r 


(3) 


where  v  is  the  mean  flow  speed  within  the  mixed  layer  and  CQW  is  the  sidewall  drag 
coefficient.  Here,  h  is  the  depth  of  the  assumed  well  mixed  layer,  which  may  be  con¬ 
siderably  smaller  than  the  maximum  depth  reached  by  the  mixed  layer  locally  at  any 
given  time.  Although  (3)  is  most  easily  derived  using  the  Boussinesq  approximation 
and  incompressibility  condition,  it  can  be  shown  to  hold  with  equal  accuracy  even  when 
the  complete  continuity  equation  is  employed.  Introducing  the  Rv  notation: 


Rv  =  bh/(av) 


(4) 


where  av  is  the  velocity  jump  across  the  edge  of  the  mixed  layer  in  the  vicinity 
of  z*h,  Price  transformed  (3)  into 
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The  implication  is  that  w  /u  should  be  two  times  larger  with  the  2LS  con- 

Q  * 

figuration  than  with  the  SOL,  for  a  given  Ry,  in  close  agreement  with  the  obser¬ 
vations.  Thompson  (1979)  obtained  the  same  result  independently,  and  had  already 
used  (7)  in  the  theory  of  Pollard  et  al.  (1973);  both  stressed  the  importance  of 
sidewall  drag  in  limiting  the  flow  speeds  induced  by  the  rotating  screen  and  in 
causing  we/u^  to  decrease  more  rapidly  than  RT  as  h/ar  became  appreciable.  The 
value  deduced  for  Ry  was  in  the  vicinity  of  0.6. 

The  present  study  was  motivated  by  the  desire  to  determine  in  future  studies 
how  the  entrainment  rate  is  modified  by  the  presence  of  a  destabilizing  buoyancy 
flux,  F^,  at  the  surface  when  a  surface  shear  stress  and  a  velocity  jump  across  z=h 
are  also  present.  The  limiting  case  would  be  free  convection  in  the  absence  of  u, 
and  Av,  for  which  Rv  =•».  An  alternate  configuration  also  producing  Rv  *«• involves 
the  use  of  an  oscillating  grid  rather  than  a  rotating  screen  to  generate  the  tur¬ 
bulent  boundary  layer.  If  Ry  can  increase  from  0.6  to  m  by  substitution  of  a 
different  mechanism  of  turbulent  mixing,  it  could  lie  anywhere  in  between  these 
two  limits  if  all  three  mechanisms  were  present  to  maintain  the  turbulence.  It 
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may  therefore  be  questioned  if  R v  is  sufficiently  constant,  during  entrainment 
that  is  driven  by  2  of  these  3  mechanisms  (u#  and  Av)f  for  (8)  to  be  valid.  To 
check  this  possibility  requires  entrainment  measurements  in  which  av  is  measured. 

Assumption  (7)  has  also  been  challenged  by  Kantha  (1978)  who  utilized  the 
model  w  /u  m  IL,  R„  ,  with  R„  being  allowed  to  vary. 

Another  reason  for  desiring  to  check  (8)  is  that  in  geophysical  situations 
the  term  will  be  absent,  yet  a  similar  equation  based  on  assumption  (7) 

(e.g.  see  Pollard  et  al. )  could  be  derived  in  which  other  forces  would  appear 
that  could  tend  to  alter  Av;  i.e.,  the  horizontal  pressure  gradient  and  coriolis 
forces.  Thus  (8)  would  then  imply  that  large-scale  forces, which  may  cause  Av 
to  change  slowly  with  time,  would  directly  affect  the  entrainment  rate.  This  con¬ 
trasts  with  the  usual  viewpoint  that  only  parameters  appearing  in  the  turbulence 
kinetic  energy  (TKE)  equation,  or  closely  allied  parameters,  affect  the  entrainment 
rate  (e.g.,  see  Price  et  al . ,  1978).  In  particular,  the  TKE  equation  does  not  con¬ 
tain  the  coriolis  parameter  but  does  contain  the  buoyancy  effect?  (8)  does  not 
contain  the  latter  except  in  self-cancelling  form.  Hence,  there  are  good  reasons 
for  questioning  the  general  validity  of  the  Rv=const  assumption  which  led  to  (8) . 

In  the  present  study  we  therefore  examine  assumption  (7)  by  means  of  direct 
measurements  of  av  and  other  relevant  quantities,  and  attempt  to  fit  measured 
values  of  we/u#  by  an  entrainment  relation  which  is  qualitatively  consistent  with 
the  TKE  equation. 


2.  Experimental  apparatus  and  procedures 

Our  primary  annulus  has  inner  and  outer  diameters  of  0.62  and  1.18m,  respec¬ 
tively,  giving  a  gap  Ar  of  0.183m.  It  is  about  0.8  of  the  size  of  the  one  KP  and  KPA 
used,  and  1.7  the  size  of  the  one  used  by  Kantha  (1978).  A  second,  inner  annulus 
was  also  occasionally  utilized;  its  outer  wall  is  the  same  structure  as  the  inner 
wall  of  the  primary  annulus,  and  its  gap  is  At  =  0.099m.  The  walls  are  very  smooth 
and  their  separation  constant  to  within  *  2mm.  A  photograph  of  the  apparatus  is 
shown  in  Fig.  1.  The  annulus  is  situated  within  a  penetrative  convection  tank 
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previously  used  in  free  convection  studies.  The  rotating  screen  is  located  at  the 
bottom  of  the  annulus,  rather  than  the  top,  because  it  was  planned  to  add  an  upward 
directed  heat  flux  through  the  screen  in  later  studies.  The  rigid  screen  is  flat 
and  true  to  ±0.7mm,  with  a  diamond-shaped  mesh  of  size  6mm  by  21mm.  Underneath 
the  screen  are  flat  insulating  sheets  of  10mm  thickness.  A  5mm  space  between  the 
screen  and  this  surface  is  taken  into  account  when  calculating  the  salt-water  mass 
budget  of  the  lower  layer;  otherwise  the  screen  is  considered  the  z=0  level. 

Two  large  plastic  windows  in  the  outer  annulus  wall  permit  illumination 
and  visual  observations  from  the  side.  The  windows  join  the  fiberglass  sidewalls 
smoothly  with  no  abutments  to  the  flow.  A  horizontally  spread  laser  beam  could 
be  positioned  close  to  the  outer  edge  of  the  mixed  layer  (at  2=hu)  midway  be¬ 
tween  annulus  sidewalls,  to  determine  mixed-layer  heights  as  in  Deardorff  et  al. 

(1980).  The  mixed  layer  was  usually  made  visible  by  adding  trace  amounts  of 
non-fat  milk  to  the  lower  layer.  The  height  of  the  laser  beam  was  continuously 
fed  into  a  16-channel  analog-to-digital  data  system  and  stored  on  tape  along  with 
other  signals.  By  comparison  with  some  salinity  profiles,  it  was  determined  (see 
Section  4)  that  h  =  h  =  1.25h,  where  h_  is  the  greatest  depth  to  which  any  mixed- 
layer  fluid  has  penetrated  at  any  given  time.  However,  in  experiments  using  temperatur 

h  lay  somewhat  below  h_.  The  laser  beam  could  also  be  utilized  to  determine 
L  2 

the  mean  slope  of  the  entrainment  interface,  associated  with  the  centrifugal  force. 
However,  h2~h  could  not  be  reliably  estimated  by  visual  methods. 

Numerous  plyolite  particles  (0.1  to  0.5mm  diameter)  with  specific  gravity  of 
about  1.02  were  also  added  to  the  water  so  that  their  passage  between  markers  0.2m 
apart  could  be  noted  and  entered  into  the  data  system  for  conversion  to  mixed-layer 
velocities  (see  Fig.  1;  the  vertical  lines  are  located  on  both  side  walls  of  the 
primary  annulus  so  that  parallax  error  is  avoided).  For  this  purpose  only  the 
central  portion  between  side  walls  of  the  primary  annulus  was  illuminated. 

In  a  few  experiments  0.1  mm  particles  were  also  inserted  into  the  non-turbu- 
lent  layer  so  that  the  mean  flow,  if  any,  at  and  beyond  the  outermost  edge  of  the 
mixed  layer  could  be  similarly  estimated.  Due  to  momentum  transfer  by  molecular 
viscosity,  v(h2)  turned  out  to  be  an  appreciable  fraction  of  v  when  wg  was  small, 
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so  that  a  correction  had  to  be  applied  to  derive  Av  from  v  (see  Section  3). 


In  3  experiments  particle  trajectories  could  not  always  be  viewed,  as  when 
using  the  inner  annulus.  A  slightly  heavy  float  of  height  hf=9.7  cm  was  then 
tracked  within  the  turbulent  layer  to  provide  v.  The  float  was  equipped  with  a 
horizontal  ring  of  diameter  0.9  Ar  which  caused  it  to  drift  midway  between 
side  walls.  Its  velocity  relative  to  v  was  calibrated  in  the  primary  annulus  as 
a  function  of  h2/h^,  yielding  downward  corrections  of  10-20%  in  v. 

In  2  SOL  experiments,  from  which  8  entrainment  data  points  were  derived,  a 
temperature  gradient  in  water  provided  the  density  contrast,  and  the  outer  layei 
given  a  linear  stratif ication  of  about  1  C/cm.  In  15  experiments  (see  Table  )) 
supplying  64  data  points,  the  2LS  was  employed  with  salt  (NaCl)  in  the  lower 
layer  of  water,  and  fresh  or  nearly  fresh  water  in  the  upper  layer.  It  may  be 
noted  that  this  configuration  differed  from  that  of  KPA  in  that  their  fresh¬ 
water  layer  was  the  turbulent  layer.  However,  the  difference  is  not  expected  to 
affect  the  results.  in  one  2LS  experiment  of  3  data  points,  kaolin  (clay)  was 
The  friction  velocity,  u*,  was  not  measured  directly  but  was  inferred  from 
the  momentum  balance  or  screen  drag  coefficient,  as  described  in  the  next  sectio' 
In  the  2  experiments  with  temperature  as  the  stratification  variable,  it 
was  measured  by  2  vertically  traversing  thermocouples  positioned  midway  betweer 
annulus  side  walls  and  separated  90°  in  arc.  They  were  mounted  from  rods  3.2mr 
in  diameter  which  descended  into  the  water  from  above.  The  thermocouple  heigh"' 
and  temperature  signals  were  monitored  continuously,  with  their  outputs  being 
averaged  together  to  provide  a  better  measure  of  the  mean. 

In  the  experiments  using  salt  or  kaolin,  the  density  of  the  two  layers  wa 
termined  gravimetrically  from  samples  withdrawn  at  the  beginning  and  end  of  e* 
experiment.  The  relative  density  difference,  la/|/^e  ,  was  kept  below  6%,  where 
f>,  is  the  density  of  fresh  water. 

In  2  experiments  with  salt,  profiles  of  f>  were  obtained  from  an  impedance  - 
(the  support  is  visible  in  Fig.  1)  of  physical  construction  similar  to  the  one 
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designed  by  Kantha  (1978).  One  needle-like  electrode  is  "platinized"  platinum; 
a  horizontal  section  of  the  stainless  steel  probe  body  provides  the  other  electrode 
situated  at  the  same  mean  height.  The  input  ac  voltage  is  of  frequency  5  KHz  ; 
the  output  dc  voltage  is  nearly  linear  with  density  for  1.004  <  1.06  when 

NaCl  is  the  constituent  salt.  Further  details  will  be  presented  elsewhere.  This 
probe  was  also  traversed  vertically,  with  density  signal  and  height  fed  contin¬ 
uously  into  the  data  system. 

For  filling  purposes,  the  annuli  and  surrounding  tank  were  placed  in  free 
liquid  contact  by  removal  of  a  cover  plate  from  a  vertical  slot  in  each  of  the 
walls  of  the  two  annuli.  They  were  then  partially  filled  with  salt  water  of 
constant  density  to  a  depth  of  0.05  to  0.1m;  later  fresh  or  nearly  fresh  water 
was  carefully  fed  in  from  above  using  a  floating  pipe  with  many  horizontally  point¬ 
ing  orifices.  The  total  water  depth  was  0.30m.  After  replacement  of  the  wall- 
slot  cover  plates,  the  water  in  the  two  annuli  was  no  longer  in  communication 
except  for  a  6mm  gap  at  the  bottom  of  the  lower  layer  —  5mm  below  the  screen 
and  1mm  above  (the  same  screen  was  at  the  bottom  of  both  annuli).  There  were  also 
three  small  gaps  at  the  bottom  of  the  inner  annulus  inner  wall  where  drive  wheels 
in  the  central  region  (see  Fig.  1)  propelled  the  screen.  The  gaps  were  partially 
sealed  with  weatherstripping,  but  still  allowed  some  diffusion  of  water  between 
the  primary  and  inner  annuli,  and  between  the  inner  annulus  and  central  region. 

In  order  that  this  transfer  not  affect  the  salt  mass  budgets,  in  6  of  the  18  pri¬ 
mary-annulus  experiments  the  lower  layers  of  the  inner  annulus  and  central  region 
were  periodically  stirred  as  necessary  so  that  their  values  of  h  would  approxi¬ 
mately  equal  h  of  the  primary  annulus.  Otherwise,  h  grew  more  slowly  in  the  inner 
annulus  than  in  the  primary  annulus.  During  another  5  of  the  experiments  a 
smooth  plexiglass  plate  was  attached  to  the  upper  side  of  the  screen  in  the  pri¬ 
mary  annulus,  but  not  in  the  inner  annulus.  Then  only  the  central  region  was 
mechanically  stirred,  the  mixed-layer  growth  rates  in  the  two  annuli  being 
more  nearly  equal.  In  the  results  to  be 
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presented,  no  significant  difference  in  scaled  entrainment  rates  could  be  detected 
between  the  11  experiments  in  which  these  precautions  were  taken  and  the  7  earlier 
ones  in  which  they  were  not.  It  is  concluded  that  the  rate  of  diffusion  of  denser 
water  from  inner  to  outer  annulus  through  the  gap  at  the  bottom  was  too  slow  to 
affect  wg.  Its  effect  on  R-f  and  Ry,  when  detectable,  was  however  taken  into  account. 

Except  in  the  2  SOL  experiments  with  heat,  several  (2-5)  constant  values  of 
screen  speed  were  employed  in  a  series  of  smoothly  connected  steps  in  each  experi¬ 
ment.  In  about  half  these  cases  the  steps  progressed  upwards  with  time,  and  in 
the  others,  downwards.  The  length  of  each  plateau  of  constant  screen  speed  varied 
from  120  to  1000s,  depending  upon  the  entrainment  rate,  and  an  interval  of  100s 
was  allowed  between  plateaus  for  each  new  equilibrium  to  be  achieved. 

In  the  SOL  experiments,  the  screen  speed,  vs,  was  usually  increased  according 
to  vs««  t  ,  which  yielded  approximately  constant  uA  values;  occasionally  it  was 
held  constant,  which  yielded  slowly  decreasing  ut  values  (see  KP) .  Screen  rotation 
rate  was  also  continuously  logged  into  the  data  system,  as  was  time  after  in¬ 
itiation  of  each  experiment. 

Prior  to  the  start  of  each  2LS  experiment  the  screen  was  rotated  for  a  short 
period,  thereby  sharpening  the  interface.  The  mixed-layer  velocity  was  then  per¬ 
mitted  to  die  out  before  actual  measurements  commenced.  Even  with  this  procedure, 
however,  in  some  of  our  experiments  in  which  the  first  plateau  employed  a  small 
screen  speed,  results  had  to  be  discarded  because  of  a  much  greater  apparent  en¬ 
trainment  rate  during  the  first  part  of  the  first  plateau  than  during  the  remaining 
part.  In  these  instances  the  initial  interface  had  apparently  been  left  in  a 
somewhat  diffuse  condition,  so  that  the  newly  developing  boundary  layer  at  first 
deepened  anomalously  rapidly  until  the  full  density  jump  was  encountered. 

Temperature  stratification  was  used  in  only  a  few  of  the  experiments  because 
the  relatively  rapid  molecular  diffusivity  of  heat  in  water  precluded  our  studying 
the  2LS  or  achieving  very  large  Rr  values.  The  larger  R^.  values  were  of  interest 
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to  us  when  it  became  apparent  that  the  centrifugal  force  was  causing  a  substantial 
mean  interfacial  slope,  and  that  only  by  increasing  and  R^.  could  we  minimize 
this  slope. 


3.  Analysis 

The  experimental  data  logged  onto  magnetic  tape  was  computer  processed 
and,  in  most  instances,  automatically  plotted  as  in  Fig.  2.  In  this  experiment 
(No.  10;  see  Table  1)  the  screen  speed,  v^,  was  initially  increased  to  a  rela¬ 
tively  large  value,  then  decreased  in  steps. 

Sidewall  drag.  At  the  ends  of  most  experiments,  and  as  in  Fig.  2,  vs  was 
increased  until  v  was  relatively  large;  then  v^  was  decreased  nearly  continuously 
so  that  in  the  primary  annulus  it  would  match  v  there  as  closely  as  possible. 

It  was  found  that  close  matching  in  situ  could  be  attained  through  constant  visual 
comparison  of  the  motion  of  the  particles  in  the  bulk  of  the  mixed  layer  relative 
to  that  of  the  underlying  screen,as  in  Fig.  2  for  t>  2450s.  Then  u^  was  essen¬ 
tially  zero,  and  it  was  also  noted  that  i?h/<>t  quickly  vanished  after  v  was  rapidly 

s 

reduced.  Thus,  the  sidewall  drag  coefficient,  CDW»  could  be  obtained  from  (3) , 
using  the  integrated  form 


ar_  fr  -  ^2) 

DW  ‘  ^iVW 


(9) 


where  subscripts  1  and  2  refer  to  the  beginning  and  end  of  the  period  over  which 

v  =  v . 
s 

The  average  value  obtained  from  15  such  determinations  is 

CDW  =  (3-7  ~  °’4>  X  10”3  ‘  (10> 


This  value  agrees  well 
walls  adopted  by  Price 


with  the  formula  for  turbulent  flow  between  smooth 
(1979)  ; 

CDW  =  °-04  (var/^)'174 


(ID 


where  V  is  the  kinematic  viscosity.  For  v  =  0.075  ms  \  a  mean  value  during 

2  -1 

the  u^=0  tests,  for  V  =  1.0  mm  s  ,  and  for  the  constant  in  (11)  extended  to 
0.040,  (11)  also  yields  3.7  x  10  3. 

It  is  assumed  that  during  entrainment  the  sidewall  drag  coefficient  is 
still  given  by  (11)  even  though  stronger  turbulence  associated  with  the  screen 
motion  and  with  av  is  then  superimposed  upon  the  turbulence  associated  with 
lateral  shear  at  and  near  the  side  walls.  A  reason  for  not  accepting  (11)  without 
first  a  direct  check  was  the  possibility  that  sidewall  curvature  should  alter 
the  formula.  The  boundary  layer  at  the  outer  concave  wall  is  expected  to  be  more 
turbulent  l  ian  at  the  convex  wall,  since  a  laminar  boundary  layer  at  least  would 
be  expected  to  be  unstable  to  Taylor-Gor tier  vortices  at  the  concave  wall.  How¬ 
ever,  the  most  noticeable  effect  from  visual  observations  was  a  reduction  of  tur¬ 
bulence  near  the  inner  wall,  as  noted  by  Kantha  (1978)  but  not  involving  as  ex¬ 
tensive  a  region  as  reported  for  his  smaller  annulus. 

Determination  of  uft.  With  CQW  obtained  from  (11),  (3)  was  preliminarily 
2 

solved  for  u,  using  the  experimental  observations  of  hjt)  and  v(t).  As  discussed 

earlier,  h  was  taken  to  be  0.8  of  the  height,  hL,  determined  by  positioning  the 

laser  beam  near  the  outermost  edge  of  the  turbulent  layer.  Usually,  the  sidewall 

2 

drag  term  was  the  dominant  contributor  to  uA  .  After  evaluation  of  these  data, 
the  screen  drag  coefficient 


was  calculated.  Results  for  the  primary  annulus  are  shown  in  Fig.  3  as  a  function 
of  the  Reynolds  number  for  the  mixed  layer,  hv/V.  When  the  screen  was  not  covered 
by  the  smooth  plate  (see  Section  2) ,  CDg  averaged 

CDS  =  6.0  x  10  3  (rough  screen);  (12b) 


when  it  was  covered  by  the  plate  it  averaged 
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CDS  =  3.5  x  10  (rotating  smooth  plate)  (12c) 

The  value  in  (12c)  is  in  close  agreement  with  (11)  for  h-*Ar,  as  indicated 
by  the  solid  curve  in  the  figure  and  the  circle  data  points,  although  the  scatter 
leaves  it  uncertain  if  the  expected  Reynolds  number  dependence  was  present.  Be¬ 
cause  of  this  scatter,  we  utilized  the  values  in  (12b)  and  (12c)  for  C  ,  and  cal- 

UD 
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culated  u*  from  (12a)  rather  than  from  (3)  and  (11)  for  purposes  of  evaluating  we/u 
and  Rf.  Examination  of  the  data  in  final  form,  using  u^  values  derived  from 
both  methods,  indicated  less  scatter  when  (12)  was  the  method  utilized,  though 
either  method  yielded  the  same  general  results. 

Our  values  of  u^  ranged  between  0.3  and  1.7  cm  s  ^(see  Table  1).  Their 

relative  uncertainty,  based  upon  uncertainty  in  (12a)  of  3%  for  vg  and  5%  for  7, 
is  6%.  An  absolute  uncertainty  in  the  mean  of  some  *  7%  is  also  present,  judging 
from  the  scatter  of  Fig.  3  and  including  uncertainty  in  the  CDg  calculation  caused 
by  uncertainty  in  relating  h^  to  h.  The  net  uncertainty  in  ut  is  therefore  esti¬ 
mated  to  be  *  9%. 

Entrainment  rate,  was  evaluated  graphically  from  enlarged  versions  of 

figures  like  Fig.  2,  taking  advantage  of  the  fact  that  wg  was  essentially  constant, 

for  the  2LS,  for  a  constant  screen  speed.  It  was  assumed  that  i>h/it  =  0.8^h^/^t. 

Typically,  20-50  independent  measurements  of  local  h^  values  entered  into  each 

estimate  of  w  .  The  error  in  w  due  to  sampling  is  estimated  to  have  been  *  15%, 
e  e 

+  -4  —  X 

and  not  less  than  _2  x  10  cm  s 

Mixed-layer  velocities.  Numerous  measurements  of  individual  particle  mean 
speeds,  as  in  Fig.  2,  were  graphically  averaged  over  each  entrainment  period  to 
obtain  the  v  values.  The  scatter  of  estimates  mostly  reflects  sampling  error 
rather  than  uncertainty  in  individual  particle  speeds.  Most  of  the  particles 
tracked  remained  in  the  central  half,  vertically,  of  the  mixed  layer  during 
their  timing,  but  some  ranged  closer  to  the  screen  or  closer  to  z=h.  The  scatter 
provides  a  lower  limit  to  an  overall  longitudinal  turbulence  intensity.  For 
cases  involving  kaolin  or  the  inner  annulus,  tbe  buoy  motion  provided  v. 
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Buoyancy  jump.  For  the  2LS,  the  quantity  b  =  g|Afl/^>,  was  evaluated  by 
asuming  that  b(h+5)  varied  smoothly  between  initial  and  final  measured  values, 
where  ii  =  Jj  ( h+h 2 )  =  (9/8)  h  is  a  mean  boundary-layer  depth 

and  S  *  8mm  is  the  equivalent  salt-water  depth  below  and  within  the  screen. 

In  many  of  the  experiments  this  quantity  was  essentially  constant,  as  expected 
for  the  2LS.  However,  in  the  rough-screen  experiments  in  which  appropriate 
mechanical  stirring  in  the  inner  annulus  and  central  region  was  not  provided 
(see  Section  2),  b(h+J)  for  the  primary  annulus  increased  by  up  to  20%  during  the 
course  of  an  experiment.  We  estimate  that  the  relative  uncertainty  in  our  meas¬ 
urements  of  bii  is  *5%,  with  a  greater  uncertainty,  crudely  set  at  tl0%,  attached 
to  incomplete  vertical  mixing  within  the  mixed  layer.  The  same  method  of  esti¬ 
mating  b  was  employed  whether  or  not  the  turbulent  layer  was  well  mixed  in  salinity. 

For  the  experiments  using  heat,  b  =  g#AT,  where  ot  is  the  coefficient  of  ther- 
n>al  expansion  of  the  water,  and  AT  is  the  temperature  jump  defined  by  T(h2)-T. 

Here,  T  is  that  mean  mixed-layer  temperature  which  at  any  given  time  yields  as 

much  warming  in  the  inner  half  of  the  mixed  layer,  relative  to  the  initial  tern-  ! 


perature  profile,  as  there  was  cooling  in  the  outer  half  of  the  mixed  layer. 

Velocity  jump.  Estimations  of  v(h2>  described  in  Section  2  disclosed  sub-  j 

; 

stantial  velocities  which,  when  scaled  by  v,  were  found  to  correlate  with  the  quantity  ] 
-  1/2 

w  (h/i/v)  as  in  Fig.  4.  With  small  w  the  mixed-layer  momentum  has  greater  ' 

e  e  j 

time  to  propagate  viscously  beyond  z=h  without  being  overtaken  by  h  (t) .  The  j 

relationship  found  for  Av  =  vfh^-v  is 

|Av| /v  =  0.81  +  0.057  ln^th^/v)*5  J  .  (13) 

It  is  ad  hoc  and  cannot  be  expected  to  hold  beyond  the  typical  conditions  encoun¬ 
tered  in  our  own  and  similar  experiments.  (See  the  Appendix  for  an  analysis  of 
how  viscosity  and  time  dependence  are  expected  to  have  influenced  Av  in  Exps.16,  17.) 

Eq .(13)  was  utilized  to  obtain  corrected  values  in  all  cases  for  which  only  v, 
not  Av,  was  measured.  Without  this  correction,  estimates  of  can  be  substanti- 
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ally  too  small  in  annulus  experiments. 

Owing  to  the  experimental  difficulty  in  identifying  particles  very  close 
to  z=h2  and  just  outside  of  the  mixed  layer  for  tracking  purposes,  the  relative 
error  in  Av  is  estimated  to  be  +  10%,  which  is  somewhat  greater  than  suggested 
just  by  Fig.  4  and  the  estimated  5%  uncertainty  in  v. 

As  found  by  previous  experimenters,  the  flow  several  centimeters  or  more 
beyond  the  interface  was  observed  to  be  laminar  in  appearance;  velocities  were 
extremely  weak  except  in  prolonged  experiments  with  especially  small  values. 

Combined  uncertainties.  The  net  estimated  root-mean-square  uncertainty  in 
w^/u^,  assuming  independent  error  sources,  is  *17%,  except  for  values  less  than 
4  x  10  4  for  which  the  uncertainty  can  approach  +50%.  The  net  uncertainty  in 
R^  is  similarly  estimated  to  be  +22%,  and  that  of  Ry,  t25%. 


4.  Results 

Temperature  and  density  profiles.  Temperature  profiles  from  experiments 
1  and  2  are  shown  in  Figs.  5  and  6,  respectively.  Although  the  bulk  of  the  mixed 
layer  has  a  significant  vertical  gradient  in  T,  the  height  h2  can  be  clearly  de¬ 
termined  and  h  estimated  to  be  some  20  to  25%  smaller.  Because  I was  rela¬ 
tively  small  in  these  experiments  with  temperature,  w^/u^  was  relatively  large; 
this  may  have  contributed  to  the  imperfect  vertical  mixing  inside  of  z=h. 

Density  profiles  using  the  salinity  probe  were  measured  only  during  experi¬ 
ments  16  and  17,  and  are  shown  in  Figs.  7  and  8.  In  Fig.  7  the  profile  for 
t=l  17s,  for  which  wg/u#  was  quite  large  (.04 3),  shows  such  a  huge  gradient  of 
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p  within  the  turbulent  layer  that  no  jump,  Lf ,  is  evident  at  all.  Many  of  the 
density  profiles  of  Kantha  (1978)  had  this  appearance.  The  other  profiles  of 
Fig.  7  were  accompanied  by  w  /u  values  considerably  less  than  0.03  and  exhibited 
a  much  more  uniform  appearance  within  the  bulk  of  the  turbulent  layer.  In  Fig. 

8  for  experiment  17  the  density  profiles  are  very  well  mixed  in  appearance  in¬ 
side  of  z=h,  and  were  accompanied  by  still  smaller  values  of  w  /u  (see  Table  1). 

•  * 

The  temperature  profiles  in  Fig.  5  were  accompanied  by  w^u^  values  in  the  neigh¬ 
borhood  of  0.025and  are  marginal  in  exhibiting  a  well  mixed  appearance.  It  might 
therefore  be  inferred  that  the  extent  of  mixing,  perhaps  as  measured  by  hOp/tz )/Lp, 
is  a  function  of  w^/u^  (see  also  Andre  et  al. ,  1979),  where  &p/tz  is  the  mean  den- 

sity  gradient  within  the  inner  3/4  or  so  of  the  turbulent  layer.  Also,  w  /u.*0.03 

e  * 

may  be  a  kind  of  threshold  above  which  the  mixing  is  highly  incomplete.  This  ap¬ 
proximate  value  also  tends  to  separate  Kantha's  poorly  mixed  salinity  profiles 
from  the  better  mixed  cases.  However,  due  to  the  existence  of  the  second  velo¬ 
city  scale,  Iv,  which  need  not  depend  linearly  on  the  first  scale,  u#,  this 
inference  is  uncertain. 

The  mean  mixed-layer  heights  noted  by  the  laser  (near  h2>  at  the  average 
time  of  each  temperature  or  density  profile  are  also  denoted  in  Figs.  5-8. 

Although  no  velocity  profile  measurements  were  attempted,  the  visual  appearance 
of  the  numerous  tiny  particles  within  the  turbulent  layer  was  always  one  of 
qualitative  well  mixedness  for  all  values  of  entrainment  rate. 

Interfacial  slopes.  During  several  of  the  experiments  h2  was  measured  near 
the  inner  and  outer  sidewalls  of  the  primary  annulus,  and  the  associated  slope 
over  a  distance  of  0.15m  determined.  From  the  balance  of  forces  for  the  lateral 
velocity  component  the  expected  slope  is 

s  =  v2/(gr  IAP I/p.)  (14) 

where  r  is  the  mean  radius.  A  scatter  diagram  of  individually  measured  slopes 
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versus  expression  (14)  is  presented  in  Fig.  9.  Although  the  data  for  small 
interfacial  slopes  lie  close  to  the  1:1  line,  the  observed  slopes  are  unex- 
plainably  smaller  than  the  theoretical  ones  for  values  of  the  latter  exceeding 
0.10  to  0.15.  However,  slopes  exceeding  even  0.1  or  0.2  are  relatively  large 
in  the  present  context,  and  are  undesirable  at  least  for  the  reason  that  en¬ 
trainment  must  be  occurring  horizontally  (u^ias  well  as  vertically  (wg) .  Only 
the  vertical  compo’-.snt  is  of  direct  interest  here. 

overall,  the  net  entrainment,  V  ,  is  always  considered  to  act  normal  and 
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outward  to  the  mean  plane  of  the  turbulent  fluid  interface,  as  in  Fig.  10. 

Only  if^h/)r  is  small  will  ug/we  be  correspondingly  small.  Considering  that 
other  uncertainties  in  w  /u  discussed  in  Section  3  lie  around  tl7%,  we  regard 
the  uncertainty  caused  by  horizontal  entrainment  (ue®0.1  to  0.2we>  to  be  of 
comparable  magnitude  and  not  cat  se  for  rejection  of  the  data.  It  may  be  pointed 
out  that  horizontal  entrainment  does  not  occur  freely  despite  the  lack  of  a 
gravitational  restoring  force;  there  is  a  corresponding  restoring  force,  from 
the  viewpoint  of  parcel  stability,  associated  with  either  the  lateral  pressure 
gradient  or  the  centrifugal  force  (see  Veronis,  1970). 

It  is  not  understood  why  slopes  of  similar  magnitude  were  not  reported  by 
previous  investigators.  We  expect  the  inter facial  slope  to  have  essentially 
the  same  magnitude  whether  the  rotating  mixed  layer  occupies  the  upper  or  the 
lower  portion  of  a  two-layer  system.  However,  L.  H.  Kantha  (personal  communication) 
states  that  substantial  slopes  were  observed,  especially  in  his  smaller  annulus. 

The  failure  of  the  observed  interfacial  slopes  to  obey  (14)  at  larger  slopes 
may  have  been  associated  with  a  strong  decrease  of  h^-h  in  proceeding  from  the 
inner  to  the  outer  annulus  wall,  so  that(^h2/^r) <^h/dr).  Another  possibility  is 
that  it  was  associated  with  a  centrifugal  secondary  circulation  which  is  upwards 
at  the  outer  wall.  This  circulation  was  quite  evident  during  the  start-up  of 
each  experiment  when  an  interfacial  slope  is  becoming  established,  and  also  dur- 
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ing  acceleration  of  the  rotating  screen  at  other  times.  However,  in  agreement 

with  previous  observations,  it  was  scarcely  evident  during  the  measurements 

with  constant,  or  quasi-constant,  screen  speed. 

Entrainment  versus  Ry  and  R  «  Our  main  results  on  entrainment  are  shown  in 

Fig.  11  where  w  /u.  is  plotted  against  RT,  as  is  conventional.  Also  shown  is 
'  e  * 

the  data  swath  from  KPA.  The  latter  lies  above  our  data,  the  customary  explana¬ 
tion  being  that  h/Ar  is  larger  for  our  data  (our  values  of  h/Ar  range  from  0.3f 
to  1.5  and  average  0.1,  while  those  presented  by  KPA  were  smaller  than  0.5). 

Eq. (8)  can  then  be  invoked  to  show  the  effect  of  increased  sidewall  friction  in 
decreasing  w^/u^. 

However,  our  measured  values  of  are  printed  at  the  data  points  of  Fig.  11, 

and  they  indicate  how  Rv  increases  from  values  near  1  up  to  5  or  10  or  more,  as 
Rf  increases  from  50  to  1000.  Although  this  is  not  surprising  in  view  of  the 
high  positive  correlation  between  Av  and  u. ,  it  means  that  the  R  =  const  concept 
which  leads  to  (8)  cannot  be  upheld  unless  almost  all  the  data  is  rejected  on  the 
grounds  that  h/Ar  was  too  large.  On  the  other  hand,  if  h/Ar  were  as  small  as 
0.1  or  0.2,  the  turbulence  Reynolds  number  proportional  to  hu#/*  would  be  con¬ 
siderably  less  than  500.  The  data  could  be  rejected  as  not  applying  to  the  large 
Reynolds-number  regime  in  which  the  value  of  V  is  irrelevant  except  for  wall 
effects. 

We  therefore  wish  to  introduce  an  alternative  interpretation  of  the  results. 
Suppose,  as  discussed  in  Section  1,  (8)  is  not  an  entrainment  law  in  gen¬ 

eral  even  for  a  laboratory  annulus  experiment.  Assumption  (7)  then  need  not  be 
considered  valid  in  general,  and  conditions  which  cause  Rv  to  increase  beyond 
some  expected  critical  value  near  unity  need  not  be  excluded  from  consideration. 

In  particular,  the  laboratory  entrainment  results  may  be  essentially  valid  even 

2 

when  the  left-hand  side  of  (3)  is  of  minor  importance  relative  to  u#  .  Then 

Rv/Rt  =  u*2/ (Av) 2  =  u*2/<v)2  =  2CDWh/Ar  (15) 
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indicating  that  R will  increase  along  with  R^  and  faster  when  h/Ar  is  greater. 

For  any  given  band  of  RT  values  our  data  tend  to  show  this  relationship,  there 
being  a  spread  in  Ry  values  such  that  greater  Rv  tends  to  be  associated  with 
greater  h/Ar. 

The  alternative  interpretation  comes  from  noticing  that,  for  a  given  narrow 

band  of  Rt  values,  w^u^  is  highly  correlated,  inversely,  with  Rv-  This  is  to 

be  expected  (since  Rv«<  (*v)  2)  from  consideration  of  the  TKE  equation  wherein 

Av  is  a  source  of  TKE  at  the  entrainment  interface  (see  Mahrt  and  Lenschow,  1976; 

Zeman  and  Tennekes,  1977;  Price  et  al. ,  1978)  and  contributes  towards  entrainment. 

Thus  we  prefer  to  consider  w  /u  to  be  a  function  of  both  R_,  and  R  ,  and  the 

e  *  i  v 

sloping  lines  in  Fig.  11  represent  our  attempt  at  estimating  a  simple  functional 
dependence  which  best  fits  our  2LS  data.  This  dependence  is 

w  /u  =  0.33RJ^R  -1'4  .  (16a) 

e  *  Tv 

It  should  be  emphasized  that  these  sloping  lines  are  hand-fitted  estimates  of 
the  values  actually  encountered,  and  are  not  assumptions  for  different  cri¬ 
tical  R  values.  An  alternate  form  of  (16a)  is 
v 

w  /lavl  =  0 . 33R—-1R  ~°’9  .  (16b) 

e  x  v 

In  (16a),  the  R^-  ^  dependence,  if  Rv  should  happen  to  be  constant,  comes 
entirely  from  the  arguments  of  Price  (1979)  and  Thompson  (1979),  using  (5)  and 
(7),  upon  considering  the  case  when  no  sidewall  friction  is  present.  However, 
we  consider  this  dependence  to  remain  valid  even  when  considerable  sidewall 
friction  is  present,  and  the  inverse  Rv  dependence  to  account  for  the  deviations 
from  an  R^-  dependence.  Eq.(5)  does  not  yield  entrainment  information  for  the 
opposite  case  of  Rr  constant  and  Rv  variable,  since  (cVdt)Rv  ft  0  then  and  (5) 
remains  a  momentum-budget  equation  which  can  say  nothing  about  the  entrainment 
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Eq . (16a)  closely  resembles  the  entrainment  relation  proposed  by  Kantha 
(1978)  and  mentioned  in  Sec.  1,  except  that  we  find  the  Rv  exponent  to  lie 
near  -1.4  instead  of  -0.5. 

Eqs.(16a,b)  are  based  upon  the  use  of  h  as  length  scale.  If  h  =  (9/8)h 
=  (9/10) h2  is  used  instead,  (16a)  becomes 


5/u.  =  0.47  (R,)"*5  (5T1-4 


(16c) 


where  the  overbars  here  refer  to  the  use  of  h  in  the  definitions.  The  precise 
definition  adopted  fc<  the  b  undary-layer  depth  thus  has  a  substantial  influence 


(40%  in  this  instance)  n  the  proportionality  constant. 

If  sidewall  friction  were  to  become  so  large  that  associated  lateral  velocity 


gradients  near  the  side  walls  were  an  important  factor,  along  with  Av  and  u^, 

in  generating  and  maintaining  mixed-layer  turbulence,  then  we  short)  expect  w  /u 

6  * 

to  depend  also  upon  a  sidewall  friction  velocity,  w  /u.  would  then  be  enhanced, 

e 

rather  than  damped,  for  given  values  of  RT  and  Rv>  Also,  if  h/Ar  much  exceeds 
unity  we  might  expect  deviations  from  a  dependence  of  w^u^  upon  only  RT  and  Rv 
because  of  insufficient  space  for  lateral  eddy  scales.  However,  for  h/Ar  as  huge 
as  1.5  we  see  no  dramatic  evidence  from  Fig.  11  that  these  other  factors  were 
very  important,  considering  the  data  uncertainties  discussed  in  Section  3.  In 
particular,  the  5  underlined  data  points  for  the  2LS  are  from  the  inner  annulus 
where  Ar  was  only  0.54  as  large  as  for  the  primary  annulus,  and  they  fit  in 
reasonably  well  with  the  rest  of  the  data.  However,  there  is  a  systematic  ten¬ 
dency  for  data  points  at  large  Rr  and  large  Rv  to  have  somewhat  enhanced  values 
of  w  /u^  relative  to  the  overall  fitted  functional  dependence.  This  may  reflect 
excess  turbulence  energy  generation  due  to  sidewall  friction  when  (Av)  was  too 
small  to  mask  it. 

Except  for  this  tendency,  the  data  uncertainties  discussed  in  Section  3 
appear  capable  of  explaining  the  scatter  of  the  2LS  data  points  relative  to  (16a) . 
(The  uncertainty  in  u.  spreads  the  data  points  in  a  direction  normal  to  the  R  = 
constant  lines  in  Fig.  11.)  Because  of  the  scatter,  the  exponent  of  Rv  in  (16) 
is  not  yet  well  determined,  nor  is  the  second  decimal  place  of  the  proportionality 
constant  accurately  known.  The  SOL  data  in  Fig.  11  will  be  discussed  in  Section  5 
the  3  kaolin  2LS  data  points  fit  in  well  with  the  salt  data. 


5.  Comparison  with  other  experiments 

SOL  vs  2LS  experiments.  With  the  portrayal  of  w  /u  of  Fig.  11  we  no 
longer  expect  any  distinction  in  entrainment  between  the  SOL  and  2LS  experiments 
(KP  vs  KPA)  that  is  not  accounted  for  by  a  distinction  in  Rv*  In  the  SOL  experi¬ 
ment  R^.  starts  out  very  small  and  progresses  rapidly  to  larger  values,  since 


V 


20  * 


2  2 

bh  =  (q/p0  )  ['i>p/*z\ terh  /2  increases  as  h  (t) .  In  the  2LS  experiment  with  uA 
constant,  is  constant;  Ry  starts  out  extremely  large  because  of  the  finite 
initial  h  and  requires  on  the  order  of  100s  (or  u^t/h  =  10  or  20)  before  it 
reaches  a  quasi-constant  minimum  value  (Price,  1979,  Fig.  3).  When  the  assump¬ 
tion  Rv=const.  is  not  made,  it  is  not  clear  if,  in  comparing  the  two  types  of 
initial-value  experiments  at  the  same  RT,  that  the  2LS  experiment  obeying  (16) 
will  have  a  smaller  Ry  value  and  consequently  greater  entrainment  than  the  SOL 
experiment,  as  in  KPA  versus  KP. 

To  check  this  point,  a  time  dependent  numerical  model  was  constructed 
obeying  (16c),  (11),  (13)  and  (3),  with  the  left-hand  side  of  the  latter  modi¬ 
fied  to  yield  an  entrainment  momentum  flux  of  -av  w  instead  of  vw  .  The  model 

e  e 

was  utilized  on  the  KP  case  with  initial  conditions:  h  (0)  =0 . 5cm,  (  1^oq)  |  c>p/i>z  |outer 
7.67  x  10  3cm  and  with  &p/pQ  =  \(h/pQ)  (ip/dz ) outer !  and  on  5  KPA  cases  with 
initial  conditions:  ii(0)=5.4  cm  and  initial  v/a  values  yielding  RT  =  36.2, 

70.3,  150.4,  292  and  523.  In  all  cases  u^  was  1.41  cm  s  ,  Ar  was  22.8  cm  and 

v(0)=0.  Results  are  those  of  Fig.  12,  the  2LS  data  points  being  those  occurring 

when  Ry  had  dipped  to  its  minimum  value  and  was  most  steady,  for  which  h/2kr  was 
in  the  vicinity  of  0.4.  The  functional  dependence  of  (16)  is  seen  from  Fig.  12 
to  have  reproduced,  qualitatively,  the  difference  observed  between  the  KP  and 
KPA  experiments;  the  latter  are  predicted  to  entrain  at  a  rate  about  1.7  times 
greater.  However,  the  KP  experiment  is  overpredicted  at  the  larger  Rr values. 

The  explanation  for  the  smaller  Ry  values,  and  consequently  greater  w^/u^ 

values  enjoyed  by  the  2LS,  is  threefold  and  no  longer  as  simple  as  when  Ry  is 

assumed  constant.  At  an  early  stage,  the  explanation  comes  from  noting  that 

the  normalized  time  (u*t/h)  at  which  (rt>sol  reaches  (rt>2ls  is  less  than  the 

time  at  which  (R  )  has  dipped  to  its  minimum  value.  At  this  stage,  when 
V  2LS 
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sidewall  friction  may  be  ignored,  and  assuming  -&v  =  v,  (3)  integrates  to 

v/u*  =  u*t/h  •  (I7) 

Thus  v/ut  was  greater  for  the  2LS,  at  the  same  RT,  because  mixed-layer  accelera¬ 
tion  occurred  over  a  longer  time  period  for  the  2LS.  It  follows  from  definitions 
(2)  and  (4)  that  (Rv) 2LS  <(Rv)g0L  and  hence,  from  (16a),  that  (we)2I/S>  ^e’soi/ 

In  the  cases  treated  by  Fig.  12,  (Rr)goL  reacbed  36.2,  70.3  and  150.7  when 
(u^t/h) S0L  was  6.7,  9.9  and  15.5,  respectively.  The  corresponding  dimensionless 

times  at  which  (R  )  was  most  steady  were  greater:  8.0,  11.8,  and  19.0.  The 
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ratio  between  the  mixed-layer  velocities  in  the  two  cases  becomes  nearly  cubed, 

through  (4)  and  (16c),  in  its  effect  upon  the  entrainment  ratio. 

At  a  later  stage,  controlled  by  wall  friction,  (v/u.J„T_  >  (v/ui)„„,  when 

x  ^Lb  *  SOL 

(R-[-) sol  =  (R-r)2LS  because  (h//Ar)$0L  had  become  appreciably  larger  than  (h/ar)2Lg. 

Sidewall  friction  was  then  stronger  for  the  SOL  case.  When  (£-),,_..  =  292  in 

Fig.  12  and  (R  )  had  reached  its  broad  minimum  value,  h  was  12.5cm  while 
V  2.  Lb  SOL 

h2LS  was  only  7.7  cm.  This  occurred  for  (u*t/h) SOL=27- 1  and  (u^t/h) 2LS=29 ■ 3 • 
This  second  explanation  takes  over  well  before  the  first  one  no  longer  holds. 

In  both  cases,  since  (w  )  >  (w  )  ,  |av|/v  is  closer  to  unity  for  the 

2LS  because  of  the  viscous  effect,  parameterized  by  (13).  This  effect  enhances 
^Rv^ SOL^^v^  2LS  and  amPlifies  somewhat  the  entrainment  rate  advantage  of  the 
2LS  in  the  laboratory. 

With  representation  (16)  there  is  thus  no  particular  relevance  to  a  finding 
of  greater  we/u#  for  the  2LS  than  for  the  SOL,  unless  both  Rr  and  R^  are  matched 
in  the  two  cases. 

Salt  versus  temperature.  Bearing  the  preceeding  discussion  in  mind,  the  tem¬ 
perature  data  points  in  Fig.  11  lie  anomalously  above  the  salt  points,  considering 
their  larger  Ry  values.  If,  for  Rr  near  100  the  temperature  data  obeyed  (16a), 

their  w  /u.  values  would  be  roughly  a  factor  of  two  smaller  than  observed, 
e  * 
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However,  similar 

behavior  was  observed  in  the  oscillating-gr id  entrainment  experiments  of  Turner 
(1968)  and  of  Wolanski  and  Brush  (1975).  In  those  experiments  an  oscillating- 
grid  velocity  scale  replaces  u.,  and  R  is  infinite  and  irrelevant.  The  entrain- 

"  V 

ment  differential  became  small  or  vanished  as  their  equivalent  Rr  variable  be- 
came  small.  Wolanski  and  Brush  looked  for  a  turbulence  Reynolds  number  dependence, 
but  found  none.  According  to  Crapper  and  Linden  (1974)  this  may  have  been  a  Peclet 
number  (Pe)  effect,  with  Pe  being  sufficiently  large  when  using  salt  that  its 
value  was  immaterial,  but  not  when  using  heat.  Here,  Pe  will  be  defined  by 
hut/D,  where  D  is  the  (molecular)  kinematic  diffusivity  in  water.  Our  experi¬ 
ment  with  kaolin  supports  this  conclusion,  since  the  Brownian -met ion  value  of 
D  for  kaolin  (with  particle  diameters  observed  to  be  near  l/im)  is  several  orders 
of  magnitude  smaller  than  D  for  salt.  In  preparing  the  kaolin  suspension  a  small 
amount  (0.5%)  of  a  def locculating  agent  was  added  in  order  to  minimize  problems 
with  the  kaolin  settling  out.  This  agent  was  comprised  of  equal  parts  by  weight 
of  sodium  silicate  solution  and  soda  ash.  However,  even  if  this  agent  had  a  mole¬ 
cular  diffusivity  of  magnitude  near  that  of  salt,  its  molecular  diffusion  across 
the  local  entrainment  interface  could  have  had  no  conceivable  effect  upon  the  net 
density  contrast.  In  contradiction,  Wolanski  &  Brush  did  obtain  very  significant  fur¬ 
ther  decreases  in  entrainment  rate  upon  utilizing  fluid  suspension^  one  of  these 
also  being  kaolin.  McDougall  (1979)  has  conjectured  that  this  result  may  have 
been  some  effect  of  a  non-Newtonian  viscosity  of  the  suspension.  However,  we 
investigated  the  viscosity  of  our  kaolin  suspension,  using  a  falling-sphere  vis¬ 
cometer,  and  found  only  the  expected  gradual  monotonic  increase  of  apparent  viscosity  j 
as  the  density  of  the  suspension  increased  from  1.00  to  1.05.  The  results  of 

\ 

Wolanski  &  Brush  on  suspensions  therefore  remain  unexplained  in  our  opinion.  J 

The  present  results  thus  suggest  the  possibility  that  the  molecular  diffusion 
coefficient  somehow  causes  enhanced  entrainment  for  RT  values  exceeding  about  50 
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(or  Rv  exceeding  about  1.3)  when  Pe  is  of  order  10,000  or  less,  and  not  when  Pe 
is  100  or  more  times  greater. 

Single  velocity-scale  experiments.  Entrainment  relation  (16b)  may  prove 

useful  in  helping  interpret  entrainment  exoeriments  designed  to  have  only  one 

velocity  scale,  lv.  If  the  present  2LS  data  are  plotted  as  w  /|Avl  versus  R 

e  v 

and  the  ut  velocity  scale  is  ignored,  they  appear  as  in  Fig.  13.  The  inner- 
annulus  data  with  larger  R^  values,  which  were  associated  with  large  h/Ar  values, 
now  appear  anomalous.  In  those  cases,  ut/  |av|  was  substantially  larger  for  the 
inner  annulus,  indicating  the  need  for  the  u^  velocity  scale  in  addition  to  |Avl. 

In  the  pioneering  experiment  of  Ellison  &  Turner  (1959)  a  layer  of  (turbu¬ 
lent)  fresh  water  flowed  over  a  weir  and  then  over  a  stagnant  saline  solution 
(along  x) ,  the  two  layeis  havinq  a  velocity  differential  Av.  However,  there 
was  a  second  velocity  scale  in  the  upper  layer  associated  with  that  entire  layer 
reraining  turbulent  after  flowing  over  the  weir  while  the  saline  layer  was  non- 
turbulent.  This  velocity  scale  was  neglected.  Since,  moreover,  it  was  x-depend- 
ent  and  not  associated  with  surface  shear,  results  of  their  experiment  can 
scarcely  be  compared  with  those  of  annulus  experiments.  Nevertheless,  the  lower 
range  of  their  data  is  shown  in  Fig.  13  to  indicate  that  theirs  was  a  small-R^ 
experiment. 

The  experiment  of  Lofquist  (1960)  comes  close  to  being  a  rectilinear  ver¬ 
sion  of  the  annulus  experiment.  A  layer  of  salty  water  was  caused  to  flow  under 
a  stagnant,  neutral  fresh-water  layer,  and  the  entrainment  rate  of  the  former 
into  the  latter  was  deduced  and  interpreted, again  only  using  the  Av  velocity 
scale.  The  second  unused  velocity  scale  was  probably  associated  mainly  with 
bottom  and  sidewall  friction.  His  mean  data  curve  is  shown  in  Fig.  13.  It  lies 
close  to  the  presentdata,  which  is  not  surprising  since  his  flow  velocities  were 
of  comparable  magnitude  (2  -  13  cm  s  ^)  and  similar  ut  values  must  have  occurred. 
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The  small  value  of  the  proportionality  constant,  c,  incurred  in  such  an 
experiment  when  analyzed  as  we/|Av|  =  cRv  n  can  therefore  be  explained  by  n=-0.9 
and  c  representing  0.3Rr  1  as  in  (16b),  with  RT  ranging  between  200  and  1200. 

In  the  experiment  of  Moore  and  Long  (1971)  jets  with  compensating 

suction  at  the  bottom  and  top  plates  of  a  toroid  propelled  two  turbulent 

layers  of  contrasting  density  in  opposite  directions.  Their  mean  data  curve 

is  also  shown  in  Fig.  13  after  adjusting  to  contain  the  length  scale  h.  It 

is  not  clear  whether  the  most  important  second  velocity  scale  here  was  associated 

with  the  3ets  emitted  from  the  slots  in  the  bottom  and  top  surfaces,  or  with 

bottom  friction  and/or  sidewall  friction.  If  the  appropriate  second  velocity 

scale  is  taken  to  be  proportional  to  4v  ,  one  would  expect  from  (16b)  an  ap-- 
-1.9 

proximate  R^  '  dependence  for  their  data  in  Fig.  13.  The  much  less  steep 
observed  dependence,  and  the  smaller  dimensionless  entrainment  at  the  smaller  Rv 
values, may  have  been  caused  by  presence  of  a  central  laminar  region  for  R^<  3 
or  6  (see  their  Fig.  4  and  discussion  on  p.644).  The  scaling  quantities  (Av( 
and  14 pi  then  each  become  a  factor  of  2  or  more  smaller  which  would  displace 
the  curve  upwards  and  to  the  right.  This  interpretation  presupposes  that  the 
appropriate  values  of  4v  and  4^  represent  property  differences  i>.  tv  ten  well-mixed 
values  and  non-turbulent  outer-edge  values. 


6.  Summary  of  results 

It  is  found  that  in  annulus  mixed-layer  experiments  where  turbulence  and 
entrainment  are  driven  by  a  rotating  screen,  u#  can  be  determined  satisfactorily 
from  the  momentum  budget  provided  mixed-layer  velocities  are  measured.  A  check 
of  this  conclusion  came  from  replacing  the  rough  screen  with  a  smooth  surface 
and  finding  that  the  calculated  screen  drag  coefficient  was  reduced  to  the  value 
appropriate  for  a  smooth  plate.  Another  check  is  that  the  entrainment  relation 
from  the  present  study  can  satisfactorily  reproduce  the  result  of  the  KP  and 
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KPA  studies. 

The  entrainment  interface  is  found  to  slope  outwards  with  radius,  with  a 
slope  as  large  as  0.1  to  0.2  when  treating  the  smaller  Rv  values  within  our  an¬ 
nulus.  Large  slopes  are  undesirable  because  horizontal  entrainment 

is  then  a  substantial  fraction  of  vertical  entrainment. 

Judging  from  temperature  and  density  profiles,  the  turbulent  layer  is  found 

to  be  quite  well  mixed  in  appearance  provided  w  /u.  does  not  exceed  about  0.03 

e  * 

(or  w  /idvl  does  not  exceed  about  0.003).  The  thickness  of  the  transition  layer 
e 

across  which  the  density  changes,  in  the  mean,  from  the  well-mixed  value  to  the 
outer-layer  value  is  found  to  be  some  25%  of  the  well-mixed  depth,  h. 

The  inverse  Froude  number,  R^,  is  found  to  vary  substantially  in  the  ex¬ 
periments.  Only  in  the  2LS  experiment  does  it  dip  to  a  quasi-constant  minimum 
value  before  rising,  but  the  minimum  value  reached  tends  to  increase  with  in¬ 
creasing  Rf  in  different  experiments. 

A  functional  dependence  of  scaled  entrainment  on  both  RT  and  R^  is  obtained 
in  (16a):  w^/u*  =  0.33Rf  **  Ry  1*4.  The  RT  **  dependence  comes  from  the  studies  of 
Price  (1979)  and  Thompson  (1979),  using  the  additional  argument  that  no  explicit 
dependence  on  wall  friction  should  occur  unless  such  friction  generates  appreci- 
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able  turbulence  kinetic  energy.  The  Rv  ’  dependence  comes  from  the  present 
2LS  data.  Together,  the  dependence  indicates  that  both  u  and  &v  promote  w  , 
with  the  Av  effect  being  the  stronger. 

With  this  representation  the  relatively  smaller  entrainment  rate  in  the  KP 
(SOL)  experiment  than  in  the  KPA  (2LS)  experiment  is  found  in  early  stages  to 
be  related  to  the  smaller  dimensionless  time  for  the  SOL  case  at  which  results 
are  compared  in  the  two  types  of  experiments.  In  later  stages  it  is  found  to  be 
due  to  greater  relative  mixed-layer  depths  then  occurring  in  the  SOL  case  which 
increases  Rv  through  sidewall  damping  of  v. 
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Direct  observational  estimates  of  Jdvl  disclose  that  it  is  typically  sig¬ 
nificantly  smaller  than  v  due  to  viscous  transfer  of  momentum  into  the  non-tur- 
bulent  outer  layer.  The  effect  is  greatest  when  wg  is  least.  The  present  data 
were  corrected  for  this  effect,  which  causes  R to  be  substantially  greater 
than  otherwise  suspected. 

The  cases  studied  using  temperature  (large  D)  instead  of  salt  to  provide  the  inter 
facial  density  jump  suggest  enhanced  entrainment,  relative  to  that  predicted  by 
(16) .  This  result  might  be  the  same  phenomenon  observed  by  Turner  (1968)  and 
Wolanski  and  Brush  (1975)  in  oscillating  grid  experiments.  An  experiment  utili¬ 
zing  kaolin  (small  D)  to  provide  the  density  contrast  in  the  2LS  indicated  no  sig¬ 
nificant  difference  in  scaled  entrainment  from  the  experiments  utilizing  salt. 

Comparison  of  the  present  2LS  data  with  other  entrainment  experiments 
designed  to  have  only  one  velocity  scale,  Av  ,  indicates  strong  similarity 
when  we/|Avl  is  plotted  against  R^.  This  result  suggests  that  the  second  velo¬ 
city  scale  in  those  experiments  associated  with  maintaining  one  or  both  of  the 
layers  turbulent  was  also  important,  as  was  u^  in  the  present  experiments. 

Although  the  present  study  made  use  of  a  wide  range  of  h/Ar  values  in  order 
to  achieve  variable  sidewall  friction  and  a  consequently  wide  range  of  R^  values, 
in  all  other  respects  small  values  of  h/Ar  are  desirable,  along  with  small  values 
of  Ar/r  and  very  small  values  of  h/r . 
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Appendix.  Viscous  propagation  of  momentum  beyond  the  interface. 

In  order  to  explore  the  plausibility  of  the  results  of  Fig.  4  and  parameteri¬ 
zation  (13),  the  Navier-Stokes  equation  was  solved  numerically  in  the  region 
h  <  z  { H,  where  H  is  the  total  water  depth.  The  coordinate  system  was  moved  ver¬ 
tically  at  the  speed  wg  whereupon  the  tangential  momentum  equation,  averaged 
laterally  (radially),  becomes 

a>v  /at  =  w  2>v  /t>z  -  ( i/~iz)  v'w'  +  j()2v  /bz2  -  2C  v  2/dr  (Al) 

c  e  c  c  dw  c 

where  v^  is  the  velocity  component  relative  to  the  moving  coordinate  system 
based  at  z=h.  The  height  h was  assumed  given  by  1.25h,  and  (^/i’zjv'w* 
by 


(a>/}z)v'w'  =  w  A.v/h)  (-5  +  4z/h) 

(h^  zg  h2) 

(A2) 

0/>z)v'w'  =  0 

(z>  h2) 

(A3) 

where  Av  =  v  ( Ii  )  -  v  (h)  . 

c  2  c 

(A2)  interpolates  the  turbulent  frictional  force  linearly  between  its  known 
value  at  the  top  of  an  ideally  mixed  layer  and  zero  value  at  h^. 

The  wall-drag  term  (last  term)  in  (Al)  was  calculated  both  using  the  tur¬ 
bulent  expression  for  CDW  of  (11)  and  the  laminar  expression  appropriate  for  a 
parabolic  lateral  profile,  the  larger  of  the  two  drag  forces  being  used. 

At  z=h (t) ,  vc  was  taken  to  be  0.9  of  the  well-mixed  speed,  v;  at  z=H  it 

was  assumed  thaOv  /Jz  =  0.  Observed  initial  heights  were  used  for  h(0), 
c 

along  with  the  observed  values  of  wg,  v  and  u*  as  step  functions  of  time  from 
Exps.  16  and  17.  The  leapfrog  numerical  scheme  was  employed,  with  the  damping 
terms  lagged  one  time  step  (At  =  2  to  5s)  and  with  a  vertical  increment  of  0.5  cm 
The  two  numerical  integrations  were  continued  over  the  lifetime  of  the  two  exper¬ 
iments  (see  Table  1)  . 

—  J5 

Results  are  shown  in  Fig.  Al  where  lAvl /v  is  plotted  versus  we(h//v)  on 

The  results  suggest  that  molecular  viscosity  was  responsible  for 


the  abscissa. 
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a  time-dependent  upper  limit  to  (av|/v  that  can  be  significantly  less  than  unity. 

However,  comparison  with  (13)  suggests  that  this  upper  limit  was  only  reached 

during  those  portions  of  experiments  having  small  values  of  u^,iavI  and  w  .  Some 

e 

mechanism  other  than  molecular  viscosity  may  therefore  have  been  important  at 
other  times.  This  mechanism  may  have  been  wave  momentum  transport  associated 
with  external  gravity  waves  centered  near  z=h^t  and/or  internal  gravity  waves 
in  the  SOL  experiments.  The  wave  amplitudes  must  increase  with  increased  values 
of  lAvl  and  u4.  The  numerical  results  thus  do  lend  some  plausibility  to  (13). 
Although  the  abscissa  stands  in  need  of  modification  that  would  include  gravity- 
wave  effects,  (13)  is  used  here  in  estimating  Av  in  the  annulus,  when  v(h^)  was 
not  measured,  because  of  the  good  data-point  correlation  of  Fig.  4. 


v 
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V‘ 


21 


TABLE  1.  Experimental  data.  Values  represent  averages  over  the  indicated  time 
period.  No  comments  in  last  column  refer  to  the  2LS  using  salt  in  the  primary 
annulus  above  the  rough  screen,  with  Av  obtained  from  v  using  Eq.(13).  Comments: 
L ll  the  SOL  using  temperature:  [2]  inner  annulus:  [3]  smooth  plate  on  top  of 
screen  in  primary  annulus  only;  [4]  4V  obtained  from  direct  velocity  measurements 
on  both  sides  of  interface;  [ 5^  the  2LS  using  kaolin. 


Exp. 

No. 

Time 

Begin 

(s) 

End 

v  (cm 

V1) 

v  (cm 
s-1) 

IA£L 

h 

(cm) 

u .  (cm 
s'1! 

IAv|(cm 

s'1) 

Rt 

R 

V 

w 

_ 6.  Comments 

u* 

la 

200 

400 

11.3 

4.8 

0.0019 

14.6 

0.50 

3.8 

109 

1.9 

0.030 

W 

b 

400 

600 

11.3 

5.2 

0.0023 

16.4 

0.47 

3.9 

167 

2.4 

0.016 

111 

c 

800 

950 

13.5 

6.4 

0.0028 

20.8 

0.55 

5.0 

189 

2.3 

0.024 

cn 

2a 

150 

250 

9.0 

2.4 

0.0006 

8.8 

0.51 

2.0 

20 

1.3 

0.066 

[i] 

b 

250 

400 

11.4 

3.8 

0.0013 

12.4 

0.59 

3.1 

45 

1.6 

0.043 

[i] 

c 

400 

550 

11.8 

4.8 

0.0017 

15.1 

0.54 

3.7 

86 

1.8 

0.019 

CU 

d 

550 

725 

12.5 

5.5 

0.0020 

17.2 

0.54 

4.4 

116 

1.7 

0.027 

ru 

e 

725 

920 

13.8 

6.0 

0.0023 

19.7 

0.60 

4.7 

123 

2.0 

0.021 

in 

3a 

300 

700 

19.2 

11.3 

0.036 

8.8 

0.61 

6.6 

834 

7.1 

0.0013 

b 

700 

1150 

19.3 

11.8 

0.035 

9.0 

0.58 

6.4 

918 

7.5 

0.0006 

4a 

600 

1100 

19.9 

11.9 

0.057 

6.5 

0.62 

6.8 

944 

7.8 

0.0011 

b 

1200 

1700 

25.4 

14.3 

0.054 

7.1 

0.86 

8.8 

508 

4.9 

0.0019 

c 

1800 

2300 

30.5 

16.5 

0.050 

8.3 

1.08 

10.9 

349 

3.4 

0.0032 

d 

2400 

2800 

35.5 

17.8 

0.040 

11.2 

1.47 

12.3 

208 

3.0 

0.0037 

5a 

100 

400 

35.9 

19.0 

0.051 

8.2 

1.31 

13.6 

239 

2.2 

0.0080 

b 

500 

900 

30.4 

17.3 

0.036 

12.2 

1.01 

11-6 

422 

3.2 

0.0034 

c 

1000 

1500 

23.6 

14.4 

0.034 

13.4 

0.87 

9.0 

590 

5.5 

0.0017 

d 

1600 

2100 

20.6 

12.0 

0.033 

14.1 

0.67 

6.2 

1016 

12 

0.0003 

6a 

1000 

1500 

19.9 

12.2 

0.033 

7.4 

0.60 

7.4 

665 

4.4 

0.0022 

b 

1600 

1900 

25.2 

14.1 

0.030 

8.8 

0.86 

9.6 

350 

2.8 

0.0056 

c 

2000 

2400 

29.8 

16.1 

0.025 

11.9 

1.06 

11.6 

259 

2.2 

0.0079 

7a 

100 

400 

35.3 

16.8 

0.025 

10.7 

1.43 

12.7 

128 

1.6 

0.012 

b 

500 

800 

30.1 

15.3 

0.018 

15.4 

1.15 

11.0 

205 

2.2 

0.0064 

c 

900 

1300 

25.5 

12.7 

0.016 

17.3 

0.99 

8.5 

277 

3.8 

0.0026 

d 

1400 

1900 

20.5 

10.5 

0.016 

18.2 

0.77 

6.5 

481 

6.8 

0.0013 

8a 

600 

1100 

20.1 

11.9 

0.029 

7.6 

0.64 

7.3 

527 

4.1 

0.0023 

b 

1200 

1700 

24.9 

13.8 

0.024 

9.4 

0.86 

9.3 

299 

2.6 

0.0043 

c 

1800 

2070 

29.7 

15.2 

0.019 

12.6 

1.12 

11.2 

187 

1.9 

0.0095 

d 

2200 

2433 

34.8 

16.2 

0.0127 

19.4 

1.44 

12.8 

116 

1.5 

0.0153 

9a 

100 

500 

10.2 

6.0 

0.021 

6.4 

0.33 

3.4 

1209 

11 

0.0015 

b 

600 

1000 

15.3 

9.4 

0.017 

8.2 

0.46 

5.7 

646 

4.2 

0.0024 

c 

1100 

1400 

20.0 

11.9 

0.015 

9.2 

0.63 

8.0 

341 

2.1 

0.0057 

d 

1500 

1730 

24.7 

13.1 

0.0128 

11.1 

0.90 

9.3 

172 

1.6 

0.0074 

e 

1800 

2040 

29.4 

13.9 

0.0092 

15.7 

1.20 

10.9 

98 

1.2 

0.0179 

TABLE  1  (Cont'd.) 


Exp. 

No. 

Time  (s) 
Begin  End 

vc  {cm 
s-1) 

v  (cm 
s'1) 

IAP1 

h 

(cm) 

u*  (cm 
s"1) 

|Avt(cm 

s-1) 

*T 

R 

V 

w 

_ S.  Comments 

u* 

10a 

100 

250 

29.7 

12.8 

0.016 

9.7 

1.31 

10.2 

89 

1.5 

0.026 

b 

350 

550 

25.0 

11.8 

0.011 

14.7 

1.02 

8.9 

152 

2.0 

0.012 

c 

650 

985 

19.9 

10.2 

0.0096 

16.5 

0.75 

6.7 

276 

3.5 

0.0026 

d 

1000 

1550 

15.0 

8.3 

0.0092 

17.1 

0.52 

4.7 

570 

7.0 

0.0007 

e 

1700 

2275 

9.1 

4.4 

0.0090 

17.4 

0.36 

2.4 

1184 

27 

0.0004 

11a 

350 

800 

15.0 

9.1 

0.0117 

6.3 

0.46 

5.8 

341 

2.1 

0.0042 

b 

900 

1150 

20.0 

10.6 

0.0091 

8.3 

0.73 

7.7 

139 

1.2 

0.0125 

c 

1250 

1490 

24.9 

11.2 

0.0051 

15.5 

1.06 

9.1 

69 

0.9 

0.028 

12a 

100 

540 

25.0 

15.3 

0.050 

6.7 

0.75 

9.4 

584 

3.7 

0.0023 

b 

640 

1120 

29.7 

16.3 

0.041 

8.5 

1.04 

10.6 

316 

3.0 

0.0028 

c 

1220 

1520 

35.3 

17.4 

0.032 

11.3 

1.39 

12.6 

183 

2.2 

0.0067 

d 

1620 

1850 

40.3 

18.6 

0.023 

16.3 

1.68 

14.3 

130 

1.8 

0.0106 

13a 

100 

195 

25.0 

11.9 

0.0113 

8.0 

1.01 

9.3 

87 

1.0 

0.026 

b 

260 

395 

20.6 

11.3 

0.0099 

10.0 

0.72 

8.1 

187 

1.5 

0.0099 

c 

470 

960 

15.4 

8.3 

0.0094 

11.1 

0.55 

5.2 

338 

3.8 

0.0022 

d 

1065 

1930 

10.1 

5.4 

0.0093 

11.4 

0.36 

2.8 

802 

13 

0.0004 

14a 

100 

200 

30.0 

11.4 

0.0065 

9.3 

1.10 

9.4 

49 

0.7 

0.045 

U] 

b 

300 

450 

25.3 

10.2 

0.0048 

14.1 

0.89 

7.8 

84 

1.1 

0.014 

L3l 

c 

600 

950 

20.8 

8.3 

0.0043 

16.8 

0.63 

6.0 

178 

2.0 

0.0089 

m 

d 

1100 

1550 

15.6 

6.2 

0.0041 

18.1 

0.43 

4.1 

393 

4.3 

0.0033 

13] 

e 

100 

200 

21.1 

10.3 

0.0064 

8.6 

0.84 

8.4 

76 

0.8 

0.051 

[3]  ,  [2] 

f 

300 

4  SB 

17.8 

8.0 

0.0047 

11.9 

0.76 

6.0 

95 

1.5 

0.014 

t3l  ,  [2] 

9 

600 

950 

14.7 

6.4 

0.0041 

13.7 

0.64 

4.5 

134 

2.7 

0.0061 

[31, 

[2  J 

h 

1100 

1550 

11.0 

4.6 

0.0038 

14.8 

0.50 

3.0 

220 

6.1 

0.0022 

[31,! 

[21 

15a 

100 

3070 

15.4 

8.3 

0.010 

7.2 

0.42 

5.3 

400 

2.5 

0.0023 

b 

3170 

3680 

20.3 

9.6 

0.0077 

10.2 

0.63 

6.6 

194 

1.8 

0.0065 

[31,  ! 

[4] 

c 

3780 

4035 

25.1 

11.1 

0.0059 

14.3 

0.83 

8.1 

120 

1.3 

0.0151 

[3],  [4] 

d 

4120 

4255 

30.1 

11.3 

0.0044 

20.0 

1.11 

7.9 

70 

1.4 

0.021 

[3],  [4] 

e 

100 

3070 

10.8 

4.3 

0.011 

5.4 

0.50 

2.5 

233 

9.3 

0.0011 

[3l,  [2] 

16a 

100 

150 

40.0 

15.3 

0.0147 

10.6 

1.46 

12.5 

72 

1.0 

0.043 

[3],! 

[4] 

b 

215 

265 

34.6 

14.6 

0.0111 

15.0 

1.18 

11.7 

117 

1.2 

0.019 

[31,14] 

c 

350 

430 

30.0 

13.1 

0.0102 

17.1 

1.00 

10.2 

171 

1.6 

0.016 

[3]  ,  [4] 

d 

500 

830 

25.3 

11.1 

0.0096 

18.6 

0.84 

8.1 

248 

2.7 

0.0040 

[3],  [4] 

e 

930 

2620 

20.1 

8.6 

0.0088 

20.4 

0.68 

5.1 

380 

6.8 

0.0014 

[3], 

[4] 

17a 

100 

1600 

20.8 

10.8 

0.024 

10.4 

0.59 

6.6 

703 

5.6 

0.0015 

[3],j 

[4] 

b 

1740 

2060 

25.6 

12.6 

0.022 

11.8 

0.77 

7.2 

429 

4.9 

0.0029 

13],  1 

c 

2160 

2420 

30.0 

14.1 

0.020 

13.0 

0.94 

9.3 

288 

2.9 

0.0043 

[3],! 

[43 

d 

2520 

2780 

35.2 

15.2 

0.017 

15.8 

1.18 

11.5 

189 

2.0 

0.0091 

[3],, 

f4] 

18a 

980 

1540 

28.3 

12.7 

0.019 

11.8 

0.92 

9.0 

260 

2.7 

0.0066 

[3] ,  DO 

b 

1650 

2600 

24.5 

10.9 

0.015 

15.3 

0.80 

7.3 

351 

4.2 

0.0032 

c 

2800 

5200 

19.9 

8.5 

0.013 

17.6 

0.67 

5.0 

9.0 

0.0007413]  ,  f5J 
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FIGURE  CAPTIONS 

Fig.  1.  The  primary  (a)  and  inner  (b)  annulus  in  which  the  screen  (d)  rotates. 
Other  parts  include:  (c)  central  region;  (e)  plexiglass  walls  of  outer  tank; 

(f)  vertical  lines  on  window  for  timing  particle  passages  (lines  on  inner  wall 
of  primary  annulus  not  visible);  (g)  salinity  probe  arm;  (h)  lasers.  Reflections 
are  present  at  far  left  and  right. 

Fig.  2.  Computerized  output  from  Exp.  #10  as  a  function  of  time  after  onset  of 

screen  rotation.  (•)  screen  speed  in  primary  annulus;  (A)  laser-height,  h  ,  esti- 

Li 

mate  of  mixed-layer  depth  in  primary  annulus;  (+)  mixed-layer  height  estimate  in 
inner  annulus  which  was  stirred  at  intervals;  (x)  mean  mixed-layer  speed  in  pri¬ 
mary  annulus. 

Fig.  i.  Mom*- n  turn-balance  estimates  of  the  screen  drag  coefficient,  C  as  a  func¬ 
tion  of  Reyn.  Ids  number  hv/i/.  (x)  rough  screen,  2ES  experiments;  (  +  )  rough  screen 
-SUL  experiments;  to)  smooth  screen.  Smooth-surface  drag-coefficient  curve  shown 
obeys  ( 1 1 i . 

Fig.  4.  Dir*  ct  measurements  of  Iziv)  ,'v  as  a  function  of  w  (h/Vv)*5  from  Exp. *16  (x) 
*26  (•,,  ar.d  «17  (•*■).  In  » I r*  and  #17  the  screen  speed  increased  in  steps;  in  #16 
it  decreased  in  steps. 

Fig.  Temperature  profiles  from  Exp.  #1  at  indicated  times.  Horizontal  line 

segments  denote  valuer,  of  h  at  time  of  traverse  of  thermocouple.  Traverses  were 
-  1 

upward  at  1  .  _  cm  s 

Fig.  6.  Te~:eraturc  profiler,  from  Exp.  #2.  See  Fig.  r>  caption  for  further  detail 

fig.  7.  Den.:  i ty  profiles  from  Exp.  #lb  at  indicated  times.  Horizontal  line  seg- 

- 1 

merits  denote  r.  at  the  respective  times.  Traverses  were  up>ward  at  from  1-3  cm  s 

Fig.  b.  Den.  lty  profiles  from  Exp.  #17.  See  also  Fm.  7  caption. 

Fig.  q.  The  retical  interfacial  slopes,  m  (see  (141  j  vet  sun  a  j  lei  slopes  in 

**  \  mj  \  » 

primary  ar.r  -  I  j.«;  ,  Hi  •  i  . Jr. 


v 


54- 

Fig.  10.  Victor  entrainment  diagram  indicating  net  outward  entrainment  velocity, 

V  :  vertical  entrainment  rate,  w  ;  and  horizontal  entrainment  rate,  u  . 

-  e  e  e 

Fig.  11.  Experimentally  determined  entrainment  diagram  of  wg/u*  versus  RT, 
using  the  well-mixed  depth  h  as  length  scale.  Printed  values  of  R^  are  encircled 
for  experiments  using  temperature  stratification  in  the  primary  annulus;  enclosed 
in  rectangles  for  the  2LS  using  kaolin;  underlined  for  the  inner  annulus  data  using 
salt  in  the  2LS ;  unadorned  for  the  primary  annulus  data  using  salt  in  the  2LS. 
Sloping  lines  obey  (16a).  The  KPA  data  swath  is  only  approximately  placed  since 
their  length  scale  may  have  been  h. 

Fig.  12.  Entrainment  diagram  from  numerical  simulation,  using  (16c),  of  a  parti¬ 
cular  KP  experiment  (curve);  and  for  numerical  simulation  of  5  KPA  experiments 
(dots).  The  (plus)  and  (square)  symbols  represent  the  respective  KP  and  KPA  ex¬ 
perimental  results. 

Fig.  11.  Entrainment  diagram  using  only  AvI as  velocity  scale,  for  the  2LS. 

Present  results  (dots,  except  I  -  inner- annul  us  results,  K- kaol in  results).  Com¬ 
parison  ir  made  with  results  of  Ellison  t.  Turner  ;  1 9  ra  9 )  ,  Lofquist  (1960)  and 
Moore  i  Lone  (1971)  . 

!•'  i  g .  Al.  Numerical  model  result.:,  neglecting  any  gravity -wave  momentum  ti  unsport 

atxjvo  z-h  ,  t  ■  Avt  v  versa:  w  lh  , v)  5  in  simulation  ot  Exp.  *  1 6  (solid-line 
a  t 

pat*. )  and  «!'  (dashed  -  1  me  path).  Letters  denote  results  occurring  during  desig¬ 
nated  periods  in  Table  1.  Along  either  path,  turn  progresses  in  the  direction  of 
the  arrows. 
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